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Abstract 

We investigate the effects of nonlinear couplings and external magnetic field on 
the thermal entanglement in a two-spin-qutrit system by applying the concept of 
negativity. It is found that the nonlinear couplings favor the thermal entanglement 
creating. Only when the nonlinear couplings \K\ are larger than a certain critical 
value does the entanglement exist. The dependence of the thermal entanglement in 
this system on the magnetic field and temperature is also presented. The critical 
magnetic field increases with the increasing nonlinear couplings constant \K\. And 
for a fixed nonlinear couplings constant, the critical temperature is independent of 
the magnetic field B. 
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1 Introduction 

Entanglement, first noted by Sclirodinger [1] and Einstein, Podolsky, and 
Rosen [2] , is an essential feature of quantum mechanics. In the last decades it 
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was rediscovered as a new physical resource for quantum information process- 
ing (QIP). Since the entanglement is very fragile, the problem of how to create 
stable entanglement remains a main focus of recent studies in the field of QIP. 
The thermal entanglement, which differs from the other kinds of entanglement 
by its advantages of stability and requires neither measurement nor controlled 
switching of interactions in the preparing process, is an attractive topic that 
has been extensively studied for various systems including isotropic [3,4,5,6] 
and anisotropic [7,8] Heisenberg chains, Ising model in an arbitrarily directed 
magnetic field [9], and cavity-QED [10] since the seminal works by Arnesen 
et al.[3] and Nielsen [11]. Based on the method developed in the context of 
quantum information, the relaxation of a quantum system towards the ther- 
mal equilibrium is investigated [12] and provides us an alternative mechanism 
to model the spin systems of the spin-| case for the approaching of the ther- 
mal entangled states [3,4,5,6,7]. But only spin-| cases are considered in the 
above papers. Zhang et al.[13,14] investigated the thermal entanglement in 
the two-spin-1 system with a magnetic field and gave a comparison between 
the uniform magnetic field case and the nonuniform one. There the nonlinear 
couplings is ignored for simplicity. In this paper, we will investigate the effects 
of nonlinear couplings and external magnetic field on the thermal entangle- 
ment in a two-spin-qutrit system. Thus we may better understand and make 
use of entanglement in QIP through changing the environment parameters. 



2 The model Hamiltonian and the solutions 

The development of laser cooling and trapping provides us more ways to 
control the atoms in traps. Indeed, we can manipulate the atom-atom coupling 
constants and the atom number in each lattice well with a very good accuracy. 
Our system consists of two wells in the optical lattice with one spin-1 atom in 
each well. The lattice may be formed by three orthogonal laser beam, and we 
may use an effective Hamiltonian of the Bose-Hubbard form [15]to describe the 
system. The atoms in the Mott regime make sure that each well contains only 
one atom. For finite but small hopping term t, we can expand the Hamiltonian 
into powers of t and get [16], 

H^e+ J(Si ■ S2) + K(S, ■ S2)\ (1) 

where J — — K — — ^ with t the hopping matrix elements, and 

e—J — K. Us{s — 0, 2) represents the Hubbard repulsion potential with total 
spin s, a potential Us with s = 1 is not allowed due to the identity of the bosons 
with one orbital state per well, since the term e contains no interaction, we 
can ignore it in the following discussions and it would not change the thermal 
entanglement. In this paper we will confine ourself to the case oi K < and 
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J < that is relevant to the recent experiment conducted on "^^Na atoms. Here 
we mainly investigate the effects of nonlinear couplings on the entanglement. 
So \K\ ^ \J\ is assumed, the Hamiltonian Eq.(l) becomes 



H = K{S^ ■ S^y, 



(2) 



with an external magnetic field, our system is described by 
H = K{Si-S2f + B{Su + S2,), 



(3) 



where Sa (a = x, y, z) are the spin operators, K is the nonlinear couplings 
constant and the magnetic field is assumed to be along the 2;-axes. When the 
total spin for each site Sj — 1 [j — 1,2), its components take the form: 
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In order to proceed we first of all find the eigenvalues and the corresponding 
eigenstates of the Hamiltonian which are seen to be 



H ^{K + 25)|*i), H 1*2) ^{K- 25)1*2), 
H IVI/3) = {^ + B)\^bs),H 1*4) = {j + 5)1*4), 

H 1*5) = ( J - 5)1*5), H 1*6) = ( J - 5)1*6), 



5|*7)=i^|*7),5 

where 



*^^ = ^i^l*^) 



(5) 



1^1) 
1*4) 

1*7) 



|1,1),|*2) = |-1,-1),|*3)=|1,0), 
|0,1),|*5) = |0,-1),|*6) = 1-1,0), 

-L(|-l,l)-|l,-l)), 



2 + (x/3t1)2 



(|1,-1) + (1tV3)|0,0) + |-1,1)). 



(6) 



here \a,P) {a = 1,0, —1 and (3 = 1,0, —1) are the eigenstates of SizS2z- The 
density operator at the thermal equilibrium p{T) = exp{—(3H)/Z, where Z — 
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Tr[exp(— is the partition function and f3 = l/ksT {ks is Boltzmann's 
constant being set to be unit kg = 1 hereafter for the sake of simphcity and 
T is the temperature), can be expressed in terms of the eigenstates and the 
corresponding eigenvalues as 

P = ^-'EexpM£^/]|*/)(*/l, (7) 



where Ei is the eigenvalue of the corresponding eigenstate and the partition 
function Z = exp[-f ](1 + 2cosh[^] + 2cosh[^] + 4exp[||] cosh[f ]). 



3 The negativity of the system 



Here we will give the entanglement of the system by applying the concept 
of negativity [17] which can be computed effectively for any mixed state of an 
arbitrary bipartite system. The negativity vanishes (i. e. negativity is equal 
to zero) for unentangled states. For our purpose to evaluate the negativity, 
in the following, we need to have a partially transposed density matrix p^^ 
of original density matrix p with respect to the eigenbase of any one spin 
particle ( say particle A) in our two-spin system which is found in the basis 
{\a,/3),a = 1,0,-1 and /3 = 1,0,-1} as p^^ = 
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where M± = ±3 + 3cosh[^] - ^Ssinhf^], Q = + exp[4^]), 

P = I exp[-f ] (3 - (3 + ^3) exp[4^]) and m = ^±fiK. The negativity 
as the entanglement measure[18] is defined by 



N{p) 



P 



- 1 



(9) 
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where 

Ta Y 



P 



P 



Ta 



^T-JpTAj+pTA (denotes the trace norm[19] of the density matrix 
le negativity N{p) is equivalent to the absolute value of the sum of 
the negative eigenvalues of p^^. From the definition N{p) = is the 

negative eigenvalue), one can see that the maximum value of the absolute value 
of the sum of the negative eigenvalues of p^"^ may be greater than 1 as long as 
tr{p^^) = 1. For different dimensions of density matrix, the maximum value 
is different. For example, the maximum value is 0.5 for the two qubits with 
two levels. In our case, the negativity is 0.5 in the state [^7), the negativity 
is about 0.972 in the state |\E'8 ) and 0.683 in the state l^^s), so the negativity 
value of the statistical mixture of these states is less than 1. In our plots, the 
value 0.972 means maximal entanglement. 




Fig. 1. The negativity vs. the nonlinear couplings constant K for different temper- 
ature T (from top to bottom, T equals 0.05, 0.6, 1.0), the magnetic field B = 0. 
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Fig. 2. (Color on line)The negativity vs. the nonlinear couplings constant K and 
the magnetic field B, the temperature T = 0.1. 



We perform the numerical diagonalization of the density matrix and the 
numerical results of the entanglement measure are presented in figures from 
Fig.l-Fig.3. Fig.l shows the plot of the negativity as a function of the nonlinear 
couplings constant K for different temperature T when B = 0. From the 
figure we can see that the thermal entanglement increases with the increasing 
\K\ and demonstrates a different evolvement curve for different temperature. 
When B = 0, K < and the temperature is close to absolute zero, the 



5 



1 

0.8 
0.6 
0.4 
0.2 

12 3 4 

T 

Fig. 3. The negativity vs. the temperature for different magnetic field B (from top 
to bottom, B equals 0.2, 1.0, 1.5), nonlinear couplings constant K = —3. 

state l^g") is seen to be the ground state, which is a entangled state, so the 
negativity is not equal to zero (in fact the negativity is equal to 0.972) and 
it decreases with the decreasing of \K\. Only when the nonlinear couplings 
\K\ are larger than a certain critical value does the entanglement exist. For 
T = 0.6, the critical value of \K\ is smaller than that for T = 1.0. We can 
also see that the maximum value that the negativity can arrive at gets smaller 
with the increasing temperature for a fixed \K\ and the system can hold high 
entanglement in a quite broader area of the nonlinear couplings K when the 
temperature is close to zero. 

The negativity as a function of the nonlinear couplings constant K and the 
magnetic field B is plotted in Fig.2 when T — 0.1. The magnetic field plays 
a negative role for the negativity at a fixed temperature, this can be seen 
from the figure. There should exist a competition between the magnetic field 
and the nonlinear couplings \K\. The critical magnetic field for various \K\ is 
different and it increases nearly linearly when \K\ increases. The negativity 
will arrive at zero when the magnetic field has a very large value. This is very 
easily understood since the term K{Si • 5'2)^ in the Hamiltonian can be ignored 
when B is very large, thus the system is unentangled. 

In Fig. 3, the plot of the negativity as a function of the temperature for dif- 
ferent magnetic field is given. When B is close to zero, the state j^g ) is seen 
to be the ground state in which the negativity is about 0.972 at T = 0. As the 
temperature increases, rapidly decreases due to the mixing of the excited 
states with the ground state. For a higher value of the magnetic field B (say 
B = 1.5), the state |^'2) becomes the ground state and there is no entangle- 
ment at T = 0. However we may increase the entanglement by increasing the 
temperature T in order to bring the entangled eigenstates such as I'^f) into 
mixing with the ground state. We also found that the critical temperature is 
almost the same for different external field. The change in negativity as T in- 
creases from absolute zero is due to population of excited levels. These results 
are consistent with those found in our previous paper [6, 13]. 
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4 Conclusions 



We investigate qualitatively (not quantitatively) the effects of nonlinear 
couplings and external magnetic field on the thermal entanglement in a two- 
spin-qutrit system in terms of the measure of entanglement called "negativity" . 
We find that the nonlinear couplings favor the thermal entanglement. Only 
when the nonlinear couplings \K\ are larger than a certain critical value does 
the entanglement exist. Increasing the nonlinear couplings \K\ increases the 
critical magnetic field, but the critical temperature is almost the same for 
different external magnetic field at a fixed nonlinear couplings constant. 
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